Quantum Electrical Dipole in Triangular Systems: 
a Model for Spontaneous Polarity in Metal Clusters 
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Triangular symmetric molecules with mirror symmetry perpendicular to the 3-fold axis are forbid- 
den to have a fixed electrical dipole moment. However, if the ground state is orbitally degenerate and 
lacks inversion symmetry, then a "quantum" dipole moment does exist. The system of 3 electrons 
in T>sh symmetry is our example. This system is realized in triatomic molecules like Na3. Unlike the 
fixed dipole of a molecule like water, the quantum moment does not point in a fixed direction, but 
lies in the plane of the molecule and takes quantized values ifJ,o along any direction of measurement 
in the plane. An electric field F in the plane leads to a linear Stark splitting =b/to • We introduce 
a toy model to study the effect of Jahn- Teller distortions on the quantum dipole moment. We find 
that the quantum dipole property survives when the dynamic Jahn- Teller effect is included, if the 
distortion of the molecule is small. Linear Stark splittings are suppressed in low fields by molecular 
rotation, just as the linear Stark shift of water is suppressed, but will be revealed in moderately 
large applied fields and low temperatures. Coulomb correlations also give a partial suppression. 



I. INTRODUCTION 

Moro et al. (de Heer's group Q, 0) found an unex- 
pected low temperature electrical dipole, apparently not 
classical, in nearly all clusters of N niobium atoms with 
2 < N < 150. Here we use the X3 molecule to propose 
an answer to the question, what a non-classical electrical 
dipole might be. Our model electrical dipole bears a close 
analogy to the non-classical magnetic moment of a spin 
1/2 particle. It is meaningless to ask the direction of the 
moment of an isolated electron, because different compo- 
nents of the moment operator do not commute. In a field 
B the degeneracy with respect to spin degrees of freedom 
is lifted and the ground state has energy Eq = —fig\B\ 
and is characterized by having maximal projection of the 
moment \±m = (0|/?Af |0) = (J-bB along the magnetic field. 

If the electron is not isolated but in contact with 
a heat bath, the statistical average of the moment is 
Pm = fJ,Bta.uh(/j,B\B\/kBT)B. The susceptibility has 
the Curie form x = ^b/^bT. The thermal behavior at 
weak applied fields is classical 3] , with a well-defined lin- 
ear response /2m = The operational definition of a 
dipole is a Curie response of the isolated system, and 
the moment is found from \f~xT . Classically we know 
that the moment points at some angle to the field and 
precesses. Quantum mechanics forces a more abstract 
view, which replaces the direction of the moment by the 
"quantization axis." 

Usually we think of an electrical dipole moment fx on a 
molecule in classical terms. When rigidly fixed in space, 
the moment has a value which is not quantized in atomic 
units edB- Its value is an accidental result of chemistry. 
In the frame of the molecule, the components of the elec- 
tric dipole moment operator commute, so a direction 
can be defined in this frame. A thermalized vapor of 



molecules has moments randomly oriented. In an elec- 
tric field F there is a net polarization proportional to 
F, corresponding to an average dipoleper molecule X-F, 
with the Curie form \ — (J^/SksT 3]. Quantum me- 
chanics enters in the usual picture when UbT is small 
enough that the spacing of quantized rotational levels 
E Iot = (h 2 /2I)L(L + 1) is no longer negligible. The 
ground state L = has full rotational symmetry with 
zero average dipole, and the Curie divergence is cut off, 
replaced by x -> M 2 /(^ 2 /2J)[l||. 

More interesting quantum effects, including non- 
commutation of the relevant matrix elements of the elec- 
tric dipole operator, may occur when the molecule has 
an orbitally degenerate ground state, and lacks an inver- 
sion center. We consider the simplest case, with equi- 
lateral triangular (D3h) symmetry. A fixed (classical) 
electrical dipole is inconsistent with symmetry. To 
quote Townes and Schawlow, "there can be no dipole mo- 
ment perpendicular to the axis of a symmetric top." 
This statement is also true in quantum mechanics if the 
ground state is non-degenerate. However, there is a non- 
zero electrical dipole matrix in the doubly degenerate E 
representation of T)3h- This matrix has non-commuting 
Cartesian components. The moment exhibits quantum 
behavior similar to the spin of an electron. Further quan- 
tum effects suppress this moment, similar to the way 
zero-point rotation suppresses the usual molecular mo- 
ment |3j,|^,|6|, but there are very interesting differences. 
Orbital magnetic moments also occur in this model, and 
are discussed in Appendix lAl 

Our model (Sec. [Hj is a simplified picture of an alkali 
trimer like Na3. We start with a rigid equilateral triangle 
geometry, fixed in space, representing, for example, three 
atoms. All electrons except the outer 3 valence electrons 
(one per atom) are frozen. One alternate realization is 
the "R-center" in an alkali halide crystal 0. In crys- 



talline NaF, for example, the center has three adjacent 
F~ vacancies with equilateral triangular symmetry in a 
(111) plane. Each vacancy binds one electron donated 
by neighboring Na + ions. Another alternate realization 
is a quantum well with triangular shape 8j, gated to 
hold three electrons, or three identical closely-coupled 
quantum dots in triangular geometry, gated to hold one 
electron each. 

Because we want our model to be relevant to Nb^v 
clusters, we must add some additional physics. We con- 
sider all distortions of the triangle in linear approxima- 
tion (Sec. IHIjl . These have a particular importance be- 
cause of a Jahn- Teller instability and a Longuct-Higgins- 
Berry phase that inverts the order of singly and doubly 
degenerate levels (Sec. IIVL and restores the ground-state 
degeneracy. We also consider rotations in space, but, for 
simplicity, only in two dimensions (Sec. EJ. Additional 
complications coming from higher order couplings, from 
rotations in the third dimension (these do not alter any- 
thing fundamental), and from antisymmetry of the nu- 
clear wavefunction (the Na nucleus is a spin 3/2 Fermion) 
are omitted for now. We plan to discuss some of these 
in a subsequent paper [3] which will focus on the Na3 
trimer. 



II. DIPOLE MOMENT OF UNDISTORTED 
TRIMER 

The minimal toy model has three electrons which live 
in a Hilbert space of three s orbitals \n), n = 0,1,2 
located on "atoms" which define vertices of a triangle. 
This crudely represents the outer electrons of an alkali or 
noble metal trimer X3. The subsequent paper will make 
contact with the specific case of Na3. 

We define positions of X nuclei as R n , n = 0, 1, 2. We 
write R n = £lr n , where f n = f n a + u n give positions of 
nuclei in an internal reference frame of the molecule. The 
orientation of the molecular frame with respect to the 
laboratory frame is given by an SO(3) matrix l~i which 
can be parameterized by three Euler angles. We choose 
the internal coordinates x, y, z such that the molecule lies 
in the xy plane. The atoms have small displacements u n 
relative to equilateral triangular positions r n o = ro£ n . 
The equilateral vertices are labeled by the unit vectors 

where £0 lies in the molecular x direction, and 1,2 are 
rotated by ±120° around z as shown in Fig. ^ 

The electrons are described by the tight-binding model 

2 

H c \ = - ^(t n ^ l+1 c n c n+1 + c.c.) + H int , (1) 

n=0 

where c„ is the annihilation operator of an electron on 
the s-orbital of the n-th atom. Coulomb interactions are 
mostly ignored. The effects of an on-site Hubbard inter- 
action Hi n t are considered in Appendix [5] The hopping 
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FIG. 1: (a) The unit vectors £ n which locate the undistorted 
positions in the molecular frame, (b) The conventional real 
basis function i\>i and (c) ip3 for the doubly degenerate E' 
orbitals which diagonalize H e \ when the molecule is undis- 
torted. Nodal planes are shown as dashed lines. If either of 
these states is occupied once and the other is empty, the atom 
will have a dipole moment ±/io, where /io = ero/2 and ro is 
the radial distance from the center to any atom. This dipole 
is shown as the bold vector. 

matrix elements t n ,n.+i depend on the distance between 
atoms \ f n - r n+ i\. 

The interaction with the external electric field F can 
be represented by 

H F = -F-(n f l) = -(n- 1 F)-fl=-F m -^ (2) 

where the electric dipole moment in the molecule refer- 
ence frame is given by 

2 

fi = -^er n c n c n . (3) 
n=a 

The factor F m = is the electric field F as seen in 

the molecule reference frame. 

For the rest of this section we take the molecule to 
be stationary, so that the lab frame coincides with the 
molecule frame (f2 = 1), and neglect distortions (u n = 
0). Therefore £„.„+i is a constant, to. For Na3, the hop- 
ping matrix element to is about 0.4 eV [lC| . Eq. Q is 
the N = 3 case of the symmetric A-atom ring, treated in 
Hiickel theory. The states follow from A-fold rotational 
symmetry 

l fc > = 7f E e ^ w >>- ( 4 ) 

* n 

The dimcnsionless wavenumber k — —1,0,1, defined 
modulo N = 3, labels the energy levels. The energy 
eigenvalues are €k = — 2to cos(2irk/N). The lowest state 
k = 0, of energy — 2to, is bonding (has no nodes), and 
is twice occupied. The other two states, k = ±1, are 
degenerate with energy +to and hold one electron. 

Now consider the electrical dipole operator J3J). For 
a single electron in the orbital |n), the dipole operator 
has expectation value (n\fl\n) = —ef n Q. The eigenstates 
\k) have equal charge on all atoms and therefore zero net 
dipole. However, there are off-diagonal matrix elements 
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obeying a selection rule Afc = ±1, 

(k+l\ft\k) = -no{x-iy), (5) 

where fig = ero/2, and x,y are unit vectors in the 
molecule reference frame. The matrix element (k\jl\k + 1) 
is the complex conjugate of the one shown. A special 
property of the trimers with half-filled valence shells is 
that there are part-filled degenerate states with k = ±1. 
The selection rule Afc = ±1 does not exclude the ma- 
trix element (k = — l|/2|fc = 1) because -1 is the same as 
2. Thus there is an off-diagonal dipole matrix element 
in the degenerate electronic ground state manifold. The 
—p, ■ F m coupling term, Eq. (|2J), is a 2 x 2 matrix in the 
subspace \k — — 1), \ k — 1), namely 

T^ndist = _ Mq(? . ^ (6) 

where a = a x x + a y y and a x and a v are the usual off- 
diagonal Pauli matrices. The eigenvalues split linearly 
Ae = ±/io|-F|, where it is important that the component 
F mz of F perpendicular to the plane of the molecule is 
omitted (that is, |F| is (F^ x + F^) 1 ' 2 .) 

One might expect that other ring molecules with odd 
numbers of atoms, like C5H5, would also have this effect. 
Indeed, there is a dipole in the degenerate manifold \k = 
2), \k — —2) for this case (since -2 is the same as 3 modulo 
5), but this level is empty and there is no dipole in the 
partly filled \k = 1), \ k = — 1) manifold. One might also 
be surprised that the Cartesian elements of the electrical 
dipole operator are represented by non-commuting 2x2 
matrices, since these are commuting operators in a full 
theory. In the full three-dimensional space of the three s 
states, the matrices fj, x , jjb y> /j, z commute, but in the k = 
±1 degenerate subspace relevant to low energy physics, 
they do not. 

One can visualize this dipole moment by using the real 
linear combinations if>2 = (\k = 1) — \k = — l))/v2i and 
■03 = (\k = 1) + \ k = — l))/v2, the conventional E' basis 
states with amplitudes (0, 1, — l)/y/2 and (2, — 1, — l)/y6 
on atoms (0,1,2) shown on Fig. These states have 
diagonal dipole expectation values ifj-ox and off-diagonal 
elements ^y. 

It is worth commenting on the magnetic properties 
of this model. The unpaired electron has a net spin. 
This couples weakly, through the spin-orbit interaction, 
to orbital angular momentum. The undistorted molecule 
has a large orbital magnetic moment. If the discrete 
three-membered ring were replaced by a 1-d continuum 
(in the spirit of the free-electron approximation for bulk 
Na metal), the orbital moment would be /t rb = kfiBZ 
where k = ±1 is the quantum number of the outer un- 
paired electron. The free electron approximation gives 
the right symmetry and order of magnitude, but is not 
quantitatively correct. A correct calculation will give 
Morb = kjfiBZ where we estimate 7 « 0.3 since the mo- 
ment must be approximately 2ireu>A with A the area of 
the molecule, and 2ttuj is the classical frequency of elec- 
tron rotation, approximately u> = to/VSfi- In the next 



sections we find that distortions quench the orbital mag- 
netic moment. For this reason we are not concerned with 
computing 7 which would have required additions to the 
model. Unlike the magnetic dipole moment, distortions 
reduce the electrical dipole moment only by a factor of 
2. The "Bloch-sphere" gives a convenient representation 
of the magnetic and electric dipoles in this model, and is 
discussed in Appendix 1X1 

III. STATIC DISTORTIONS AND ADIABATIC 
LONGUET-HIGGINS PHASE 

To describe small distortions, consider the first order 
expansion in u n = f n — ro£ n . The hopping matrix is ex- 
panded to first order around the undistorted atom sepa- 
ration V3ro, 

t n ,n+l = h - -y=(\f f n - T n+ i\ - V3r ) 

= h - -(U n +1 - U n ) ■ (in+l ~ In), (7) 

where g is the electron-phonon coupling. 

To completely specify three atomic coordinates R n re- 
quires nine numbers. Three give the position of the center 
of mass, and describe rigid translations of the molecule 
which separate out and do not contribute to Eq. J7J). 
Three (e.g., Euler angles) locate the spatial orientation 
of the plane of the molecule, and the rotation angle of 
the molecule in the plane. These rotations, which also 
do not contribute to Eq. J7J), will be discussed later. Fi- 
nally, three parameters, chosen as amplitudes of three 
orthogonal normal modes of distortion, give the vibra- 
tions which do appear in Eq. (JJJ. 

The normal modes of an equilateral triangle comprise 
a symmetric ( A[ , or "breathing" ) mode and a doubly de- 
generate (E') mode. Their amplitudes are conventionally 
[ill IT2]] called (Qi, Q2, Qa)- Fig. |21 introduces a second 
trio of unit vectors v n . Like the previously introduced set 
£„, these vectors lie at 120° to each other; they add to 
zero and are orthogonal to all translations and rotations 
and to the breathing mode eigenvectors £ n . 

The atomic coordinates omitting rotations and trans- 
lations are 

r n = (ro + Qi)in + Qi>„ (8) 
|„ = (cos(27rn/3),sin(2 7 rn/3),0) t (9) 
v n = (cos(a — 27rn/3), sin(a — 27rn/3), 0)^ (10) 

where the dagger indicates the transpose of the row vec- 
tor, which is a column vector. The amplitude Q = 
\JQ\ + Q\ liil an d angle a — tan _1 ((53/Q 2 ), shown 
on Fig |3 parameterize the E' distortion. We will also 
assume that Q\,Q r so that molecule is close to 
its symmetric "equilibrium" (see below) configuration 
Qi = Q = or f n0 = r i n . 

Vibrations are especially important because the or- 
bitally degenerate ground state is Jahn- Teller unstable. 
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FIG. 2: Geometry of the distorted trimer. Atoms are at co- 
ordinates f n = (ro + Qi)£n + Qv n - When Q = the molecule 
has full symmetry. When the angle a = 0, the vectors v n 
lie along the dashed lines and describe a symmetric stretch 
which makes an isosceles distortion. This is the mode called 
Q2 by Van Vleck. When a — 7r/2, the vectors v n give the 
asymmetric stretch Q3. These two patterns form the con- 
ventional basis for the doublet E' vibrations. Rather than 
(Q2, Q3), we use the parameters (Q, a) to describe the E' vi- 
brations. Note that the three vectors v n lie at 120° and add 
to zero, independent of a. Note also that when a = 120°, the 
distortion is isosceles with atom 2 at the vertex, and when 
a = 240° , the isosceles distortion has atom 1 at the vertex. 



The doublet E 1 distortions lift the degeneracy, lowering 
the energy of one of the doublet E' electron states. The 
latter is singly occupied in the ground state, thus lower- 
ing the energy of the molecule. In this section we make 
the static approximation. We assume that the trimer 
is distorted so that Qi.2,3 7^ and solve electronic part 
of the problem in the presence of this static distortion. 
First, project Hamiltonian onto the doubly degener- 
ate ground state of an undistorted trimer. This projected 
Hamiltonian gives the electron-phonon term acting in the 
subspace of (\k = —1), \k = +1)) states as 

W ep = gQil - g(Q2<Tx - Qzv y )- (11) 

The "breathing" amplitude Qi shifts the degenerate lev- 
els uniformly. It also shifts the doubly occupied \k = 0) 
level. Both these shifts will be absorbed into redefining 
the size ro to be the true radius of the molecule. The 
(Q2,Q3) vibrations split the degenerate single-particle 
energies by ±gQ where Q is the magnitude of the 
(Q2, Q3) distortion. The splitting does not depend on the 
angle a. The "lower Jahn- Teller" electron orbital (which 
is the lower eigenstate of Ea. i|ll|) . occupied with unit 
probability in adiabatic approximation) has the wave- 



function 

V> ad = ^={e-^ 2 \k = l)+e^ 2 \k = -l)) 

= sin(a/2)V'2 + cos(a/2)V> 3 - (12) 

Here the first form of f/'ad uses the basis of Eq. (0} while 
the second form uses the basis of real functions shown 
in Fig. H For a — the distortion is an acute isosceles 
triangle and the lower Jahn- Teller state is ip 3 as shown in 
Fig. ^ panel (c). For a = 180°, the distortion is obtuse, 
and the lower state is ip2- The node of the lower state 
cuts the longer bond, leaving the shorter bond nodeless. 

The effective interaction with the external electric field 
can be obtained from H6I12|) as 

n F = (^ ad |H r / idgs |v ad ) 

= — fiaFm ■ (—x cos a + y since). (13) 

The state ipad has an electrical dipole (ipudlfilip&d) = lk)P, 
where the unit vector P is {—Q2X + Q3y)/Q- A correc- 
tion of order (Q/ro) cos(3a) to the dipole will also occur 
because the values of the vectors r„ in Eq. change 
under distortion. However, effects ignored in our model 
make similar size corrections, so we ignore this. The 
dipole magnitude fio is the same as for the undistorted 
Dsh case, but it points in a fixed direction. For a = 0, the 
distortion is acute and the ground state wavefunction the 
same as Fig. ^c), with fl in the —x direction. As a in- 
creases in a counter-clockwise sense, the negative charge 
of the last electron circulates clockwise (corresponding to 
current moving counterclockwise) . The dipole fl counter- 
rotates as a rotates. We have lost the quantum nature of 
the dipole, which now orients classically. However, quan- 
tizing the distortional degrees of freedom will restore the 
quantum aspect. The orbital magnetic dipole is quenched 
in the adiabatic ground state, since the k — ±1 compo- 
nents are included in -0 a( j with equal magnitude. Quan- 
tizing distortions will restore the magnetic dipole only 
very weakly. 

Notice that the adiabatic wavefunction Eq. i|12|) 
changes sign when the molecule is taken once around 
the point Q — of "conical intersection" of the upper 
and lower Jahn- Teller energy surfaces. That is, when 
a increases by 27r, -0 a d changes sign. This adiabatic 
phase 7T picked up by -0 a d is known as the Longuet- 
Higgins phase^4| or as Berry's phase^f|- Important con- 
sequences of this for the quantized distortional states and 
for the electric dipole are worked out in the next section. 



IV. QUANTIZED VIBRATIONS ALTERED BY 
LONGUET-HIGGINS PHASE 

The full Hamiltonian of the molecule can be written as 
W = H ions + iJpot + H c \ + Hp, (14) 
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where 



h 2 JL, d 2 



Alans - 9Af E 



2M ^ dR 2 ' 



(15) 



is the kinetic energy of the Na + ions. The potential en- 
ergy is approximated by central springs, 



Hp t = ^-r(|fn-r n +i|-V^o)' 

K 
~ ~2 



3Q? + 2^ 



(16) 



We use cyclic notations n + 3 = n; H e i, Hp are Hamil- 
tonians 1)112(1 of the valence electrons and of coupling to 
an external electric field F. 

In adiabatic approximation one looks for the total wave 
function of electrons and ions in the form of a prod- 
uct ^(R^tpad- Then the effective Hamiltonian acting 
on ^j(R n ) is obtained from Eqs. (|1 1115116(1 . It acquires 
an induced "potential energy" (-0ad|H cp |V ; ad) = gQi—gQ 
and can be written in the reference frame of the molecule 
as the sum of the kinetic and potential energy of vibra- 
tional modes, 



Ha = 
He' = 



Ha 
3M 
~ 
3M 

~ 

3M 

~2 



H, 



Qi 



Qi 



gQi 

2 ^ 1 6K 

v Qi) + -rQ 2 - gQ, 



(17) 
(18) 



where M is the mass of the Na atom. The second 
form of Ha, Eq. i|17|l . results from the redefinition 
r o — g/^K —* r Q of the size of the trimer, and the corre- 
sponding redefinition Q\ + g/3K — > Q±. The breathing 
mode frequency u>a in this model is y/ K/M; in Na3, mul- 
tireference configuration interaction calculations give this 
to be about 17 meV 0], while experiment suggests 
16 meV, which implies K 1.5 eV/A 2 . 

Unlike Eq. fT7|l . Eq. fT%|) cannot be simplified by 
a shift of origin. Rather than being a simple shifted 
parabola, the two-dimensional potential surface of Eq. 
(II Hp has a cusp at the origin. This is the "conical intersec- 
tion" between the lower energy adiabatic energy surface 
and the Jahn- Teller-split higher energy surface, which 
replaces — gQ by +gQ in Eq. JTSJ. The higher surface 
does not enter in adiabatic approximation. The poten- 
tial energy of Eq. (|18fl is cylindrically symmetric around 
Q2 = Qi = Q = 0, and has the form of a "Mexican hat." 
The value Q m i n = 2g/9K minimizes this potential in the 
radial direction in (Q2, Q3)-space, and the energy is low- 
ered by the Jahn- Teller energy —g 2 /9K. The frequency 
of radial oscillation is uq — y/3K/2M, but for angular 
motion there is no restoring force, uj a = 0. The free mo- 
tion of the variable a is called "pseudorotation." When 



the coordinate a increases monotonically, the molecular 
coordinates r„ of the atoms move in counter-clockwise 
paths in the molecular x, y plane (see Fig. |2J) accompa- 
nied by no collective rotation, but by true angular mo- 
mentum, which of course must be quantized. 

Quantum chemical calculations for Na3 give a some- 
what different picture, because the distortion (~ 10%) 
is not very small, and as a result there is significant 
warping of the path of least energy on the energy sur- 
face in the Qi , Q3 plane surrounding the conical inter- 
section at the origin. Most calculations find the min- 
imum to be an obtuse isosceles triangle (a = 180° or 
equivalently a = ±60° on Fig. |2J - Saddle points, rep- 
resenting barriers for pseudorotation, occur at the acute 
positions (a = 0° or equivalently ±120°) [l^. The calcu- 
lation of von Busch et al. [l6| gives the average distortion 
Qmin = 2g/9A" sa 0.25 A. which implies that the coupling- 
constant g is about 1.7 eV/A. The same calculation gives 
the energy minimum to be 97 meV and the barrier height 
to be 25 meV, higher than the breathing phonon energy, 
and sufficient to strongly inhibit pseudorotation. In our 
model the energy minimum — g 2 /9K is about 200 meV, 
twice the quantum-chemical value, and there is no barrier 
in leading order. 

The kinetic part of Ea. (|18|) can be written as 
(3M/2)Q 2 + (/q/2)q; 2 , in separated radial and pseudoro- 
tational parts, with pseudorotational moment of inertia 
Iq = 3MQ 2 . The zero-point amplitude of radial oscilla- 
tion (H 2 /6KM) 1 / 4 is small enough compared with Qmin 
that we can approximate the pseudorotational moment 
of inertia by the constant Iq = 3M Q 2 nin - Then the pseu- 
dorotational wavefunction (in our approximation with no 
restoring force) is ^ps-rot = exp(imo;) and the energy is 
H 2 m 2 /2In. As pointed out long ago by Longuet-Higgins 
et al. [I9j, the pseudorotational quantum number m 
must take half-integer values. This is because the total 
wavefunction ip pa - I0 t(Q, oOV'ad must be a single- valued 
function of the coordinate a. We already noticed that the 
adiabatic electron wave-function Eq. H12|) changes sign 
under one pseudorotational cycle. Therefore the function 
"0ps-rot must also change sign when a is increased by 
2n. The ground state, with pseudorotational quantum 
to = ±1/2, is thus still doubly degenerate. The elec- 
tronic orbital degeneracy was lifted by the Jahn- Teller 
distortion Q — > Q m i n , but left and right pseudorotating 
states restore the degeneracy of the combined electronic 
and atomic system. This restoration of degeneracy is the 
discrete system analog of restoration of translation in- 
variance by phase solitons in one-dimensional molecular 
chains poj . 

The wavefunctions ip m of this adiabatic vibrational 
ground state are products of the pseudorotational part 
exp(zTOa) times the electronic part -0 ac i (Eq.JT^J, also 
multiplied by the ground state harmonic oscillator func- 
tions for breathing Q\ and radial oscillation Q. The 
dipole operator in this 2d manifold is obtained by pro- 
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jecting Eq. (|13fl onto 



±1/2) 



dipole (PQD) is 



n 



ad vib gs 



(19) 



where f = t x x ± T y y are Pauli matrices acting in the 
space of pseudorotational states \m — 1/2) and |m = 
— 1/2). The eigenstates ±1/2 of pseudo-angular momen- 
tum (—id a ) have left and right circulating dipolar charge 
currents. Unlike the undistorted molecule, the frequency 
of circulation is no longer to/h, but is now fi/Jg, a much 
smaller number. Therefore the orbital magnetic dipole 
moment is restored, except with a much smaller value. 
There is no diagonal expectation value of the electric 
dipole, only an off-diagonal value, just as in the undis- 
torted D 3 h case, but reduced by 2, from /i to fJ-o/2. 
Thus the half-filled outer orbital of the symmetric tri- 
angular molecule has a quantum dipole which persists 
under quantization of molecular distortions. A field F m 
applied in the plane of the molecule splits the levels lin- 
early ±fi \F m \/2. The eigenstates in the field are linear 
combinations of the left and right pseudorotators whose 
moment is distributed in the plane of the molecule with 
probability distribution cos 2 (o:f) where is the angle 
measured from the projection of the field onto the plane 
of the molecule. This cosine distribution is what causes 
the reduction of fj,Q by 2. 

In principle this is not fundamentally altered by a 
warping potential like (Vb/2)(1 + cos(3a)), because tun- 
neling along the pseudorotational minimum energy path 
preserves 3-fold symmetry and guarantees a degenerate 
ground state [2lJ ■ Pseudorotational quantum numbers m 
are still half-integer quantized, but defined modulo 3, so 
that for example m = 3/2 and m = —3/2 are the same 
and this state is singly degenerate, while m = ±5/2 is a 
doublet belonging to the m = ±1/2 representation. The 
splitting fiwpseudorot between the m = ±1/2 ground state 
and the m = 3/2 excited state gets reduced in frequency 
from uj — H/Iq by a Franck-Condon-type overlap factor 
governing the tunneling probability. Thus, warping can 
greatly reduce the orbital magnetic dipole, which is pro- 
portional to the pseudorotational splitting. Also, thermal 
mixing of pseudorotational states will smear any effect for 

ksT > fiWpseudorot- 



V. QUANTIZED ROTATIONS 

Now we must deal with real rotations of the molecule in 
space, which adds a lot of complexity, as it does even for 
a fixed dipole. A full 3-d analysis does not fundamentally 
change things from a much simpler analysis of rotations 
in 2-d, so we confine ourselves to 2-d rotations in this 
paper. Here we describe how 2-d rotations alter both a 
fixed dipole and our quantum pseudorotator. We work 
in the limit where the molecule stays in its vibrational 
ground state, but has pseudorotational excitations. 

The Hamiltonian for our pseudorotating quantum 



^2dPQD = + + y(l±cos3a) 

- H Fcos(ip + a) (20) 

Here ip is the Euler angle defining the orientation of the 
reference frame of the molecule with respect to lab frame. 
7r a and tt^, are dimensionless angular momenta which in 
quantum theory are — id a and — id^ and take half-integer 
and integer values correspondingly. The first term of 120|) 
is the kinetic energy of rotation of the trimer as a whole 
with Iq = 3M 7q . The second term is the kinetic energy of 
pseudorotation, and the third term is the simplest model 
for the warped potential energy of pseudorotation, with 
minima at the obtuse isosceles positions, and barriers of 
height Vq at the acute positions. The simplest version of 
our model has Vq = 0. The last term is the coupling to 
the external electric field. It can be obtained from l|13|) 
by substituting F m = F (cos tp, — simp). 

Eq. with V set to zero, models a quantum elec- 

tric dipole moment, free to rotate in 2d space. We now 
calculate the observable dipole moment for a statistical 
ensemble at temperature T. This is (/z) = —dF/dF, 
where T is the free energy of the dipole calculated at 
temperature T with Hamiltonian (|20|l . The temper- 
ature scales set by rotations and pseudo-rotations are 
T rot = h 2 /2I < T ps _ rot = h 2 /2I Q . For Na 3 these are 
approximately 0.1K and 5K. We assume T ^> T Iot , which 
allows us to treat spatial rotations classically, but T can 
be either above or below T ps _ rot . 

Before proceeding to the quantum dipole l|20|l we first 
summarize the well known results for the classical dipole 
moment in an external electric field. 



A. Classical dipole 



The Hamiltonian of a 2d fixed dipole is 
Ti-2d FD = ^aFcosip 



2/ 



(21) 



where Jo is the molecule's moment of inertia. This is the 
pendulum problem. The low-lying quantum eigenener- 
gies of this system, calculated numerically, are shown in 

Fig. m 

The classical treatment is standard. After evaluating 
the partition function Z = exp(— J-/T) (we use units 
k B = 1), the classical average dipole is found by differ- 
entiating T to be 



(m) = Mo 



J* dlpe^oF/T) cos v 



(22) 



which saturates to (1q at high field F and equals xF a t 
low field, where x = Mo/^T is the 2-d version of the 
Debye-Langevin polarizability [3j. 
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FIG. 3: Spectrum of low-lying energy eigenvalues of the 2- 
d fixed dipole rotator. The field energy is fioF and the 
rotational quantum of energy is T ro t = h 2 /2I. The zero 
field axis has dimensionless eigenvalues equal to m 2 for m = 
0, ±1, ±2, .... Each level except m — is doubly degenerate, 
and the degeneracy is split by the field. 



FIG. 4: Dipole moment of classical and quantum dipoles as 
a function of applied electric field, calculated at temperature 
T = 7.0Trot. The quantum rotator has a dimensionless pa- 
rameter p = TQ/T ro t which is set to 40.5 (dashed curve, rea- 
sonable for Na3) and to 200 (dotted curve). 



The results of the quantum treatment of this problem 
do not have a simple form. At high field (hqF 3> T) the 
dipole makes pendulum oscillations around the direction 
of the field. This is seen in the lower right part of Fig. [3] 
where there are evenly spaced levels at AE = hu> where 
the pendulum frequency lu — y/ JioF/Iq increases as the 
field increases. At low field and low temperature, the 
quantum ground state is rotationally symmetric with a 
dipole induced by second-order coupling to the first ro- 
tationally excited states. In this limit the polarizability 
becomes jUjj/ (^ 2 /2-0- If the field increases to values large 
compared with T/fiQ, the eigenstates are mixtures of ro- 
tational levels. These mixtures give a full dipole moment, 
saturated at the value (fi) — fiQ. We plot the dependence 
of dipole moment versus an applied electric field in Fig. 
21 In different limits we have 

„ !*L « Trot ) 

MO J- rot 

~ m o t « « T) 

* i-^-B (23) 

MO-T 1 rot 

where the subscript "CD" means classical dipole (the ro- 
tations are treated quantum mechanically). The water 
molecule is a good example. Normally we think of a 
classical fixed dipole, but the quantum eigenstates in zero 



field are rotating dipoles with no fixed average moment. 
In a typical chemical or biological environment, the en- 
vironmental field F is large, free rotations are hindered, 
and the dipole is exposed. 

B. Quantum dipole with free pseudorotations 

To emphasize the difference between a quantum dipole 
(|20[) and a fixed dipole (|21[1 . we choose the limit Vq = 0, 
i.e., a quantum dipole with free pseudorotations. Eq.J^OJ) 
with Vq = is a funny degenerate sort of double pendu- 
lum with two angular degrees of freedom. The quantum 
states at zero field F = have eigenvalues 

e(£, m) = h 2 (£ + m) 2 /2I + h 2 m 2 /2I Q (24) 

where the angular momentum quantum number £ takes 
integer values 0, ±1, . . . and the pseudomomentum quan- 
tum number m takes half-integer values ±1/2, ±3/2, . . .. 
Because Iq -C Io, the low- lying states have m = ±1/2. 
All eigenstates are 4-fold degenerate, even the ground 
state, because for either choice of |m| there are two ways 
to get any allowed value of \£ + m|. When the field is 
applied, each quadruplet level splits into two doublets. 
There are no allowed matrix elements of the field in any 
degenerate manifold of 4 states, so the initial splitting 
is quadratic in field, just as for water. If we assume 
that the temperature is low enough that only the low- 
est pseudorotational levels m = ±1/2 are excited, it is a 
good approximation to ignore the term in Eq. I|2U|) involv- 
ing Iq. Then the 2-d Hamiltonian (|20|) separates into 
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FIG. 5: Spectrum of low-lying energy eigenvalues of the 2- 
d quantum dipole rotator. The field energy is jj,oF and 
the rotational quantum of energy is h 2 /6MrQ. The zero 
field axis has dimensionless eigenvalues equal to (£ + m) 2 + 
{ro / q-min) 2 m 2 , with the zero-point pseudorotational quantum 
(fo/'i , min) 2 (l/2) 2 subtracted out. The level in red at zero F is 
the first with pseudorotational quantum number jraj > 1/2. 
Each level at zero F is quadruply degenerate, and the degen- 
eracy is reduced to double by the field. 



2x2 matrices acting on subspaces formed by \l,m = 
1/2), \l — l,m — —1/2). The ground state energy is 
r ro t [5/4- y/l + (fi F/(2T rot )) 2 }, which implies a ground 
state dipole (fi) = {^/2)[(^F /2) / V( Mo F/2) 2 + 7*J. 
For high fields and T not too low (/j, F > T > T rot ) the 
dipole is 



</'■> 



2d QD 



Mo 



2T'T ro t 



(25) 



Notice that the effect of temperature (when it is less 
than the field energy (J-qF) is smaller on the quantum 
(pseudo- + real-) rotating dipole Eq. l|25|l than on the 
fixed rotating dipole Eq. H23J1 ; there is an extra small fac- 
tor T rot /(p, F/2) which postpones thermal degradation 
of the dipole by rotation. The reason is evident in the 
difference between the spectra seen in Fig. [3] and Fig. 
|SJ At higher fields F and not too high excitation energy, 
the fixed dipole has energy levels that fan out to larger 
separations in higher fields, meaning lower slopes and 
lower dipoles in higher excited states, while the quantum 
dipole has energy levels whose slopes do not depend on 
excitation level. The fixed dipole levels correspond to 
pendulum oscillations of the dipole away from the field 
direction, whereas the quantum dipole levels correspond 
to free rotations (the energy levels are H 2 (£ + l/2) 2 /2/ ) 
of a molecule whose dipole is fixed in the direction of 
the field. The way this happens is that the dipole is not 



fixed in the axes of the molecule, but free to pseudoro- 
tate around these axes. Real rotations and pseudorota- 
tions give compensating changes in the dipole moment, 
leaving the moment fixed in space. This happens equally 
well for left and right rotations, so all levels are still dou- 
blets. The comparison between dependences of (/z) vs. 
F is presented in Fig. 0] 



VI. DISCUSSION AND CONCLUSIONS 

The main point of the present paper is the possible 
relevance to the issue of unusual polarity in metal clus- 
ters pi ■ Existing experiments on X3 molecules have 
not seen the effects described here. There are multiple 
possible reasons, including (a) not looking for quite the 
right effect, (b) need for lower temperature, (c) nuclear 
exchange symmetry, and (d) suppression by large barriers 
to pseudorotation. A subsequent paper will address some 
of these issues in the context of Na3. Here we speculate 
about the more general issues. 

Suppose we discard triangular symmetry altogether, as 
in a mixed alkali trimer like L^Na [2^. or LiNaK. In the 
symmetric case, symmetry forced the conical intersec- 
tion which preserved (by pseudorotation) the degenerate 
ground state and its quan tum dipole. However, as ex- 
plained by Wigner [T3.l23j| , points of degeneracy, or con- 
ical intersections, are generically present [24[ if there are 
enough free variables. In the last occupied molecular or- 
bital level of LiNaK, the splitting of the two sublevels can 
be tuned to zero by adjustment of two parameters, since 
the Hamiltonian matrix can be written in this manifold 
as real symmetric 2x2. There are three such parameters, 
the normal mode amplitudes (Qi, Q2, Qs)- Thus it is per- 
mitted to have a line in configuration space where the up- 
per and lower energy surfaces touch and away from which 
they separate linearly. Consider the minimum energy 
point on this line, and also the global minimum of the 
lower energy surface. Away from the global minimum, 
the energy increases quadratically with three principal 
axes defining nearby ellipsoidal energy surfaces. Now be- 
gin moving along the axis of least energy increase, and 
continue following the bottom of this valley. It is possi- 
ble that this valley will return to the global equilibrium 
after circling the line of accidental degeneracy. It is also 
possible that the path in this minimum valley lies always 
lower than the minimum energy found on the line of de- 
generacy. In this case, there is again a pseudorotation 
with an electronic-wavefunction sign change. The low- 
est energy pseudorotational state must therefore have a 
compensating phase change of it when orbiting the degen- 
erate line. Time-reversal symmetry then requires that an 
inequivalent reverse pseudorotation with opposite phase 
is degenerate. These two states should each drag a dipole 
moment in opposite directions once around the degener- 
acy line as the system moves adiabatically around the 
circuit. Thus a quantum dipole is generically permitted. 
Larger clusters of metal atoms have more complex en- 
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ergy surfaces, and both features (1) conical intersections 
and (2) lack of inversion symmetry are generically likely. 
Therefore unusual dipoles whose properties are not fully 
classical may be expected. Possibly this is what Moro et 
al. P10 observed. 

Moro et al. observed a correlation between unusual 
dipoles and metals known in bulk form to have super- 
conductivity caused by large electron-phonon coupling. 
Large electron-phonon coupling derives from a high elec- 
tronic density of states and a rapid change of single- 
particle energy eigenvalues with deformation. These are 
conditions favorable to the occurence of the dipole we 
have described. A high density of states means plentiful 
levels whose energies may accidentally cross at conical in- 
tersections. Sensitivity of eigenvalues to distortion means 
that the point of the cone of intersection can lie higher 
in energy than the surrounding landscape. Thus our pic- 
ture contains the right ingredients. Further experiment 
and more accurate microscopic modelling will be needed 
to confirm it. 

To conclude, we studied the formation of an electri- 
cal dipole moment in triangular molecules. Although, 
the classical dipole moment is forbidden for symmetric 
trimers it is formed due to the quantum degeneracy of 
the ground states. This dipole moment survives Jahn- 
Teller distortion of the molecule and keeps its "quantum 
nature" for small distortions. It is interesting to note 
that the quantum dipole moment makes a two-level sys- 
tem which can be manipulated by an external electric 
field. If it were possible to have a controlled interac- 
tion between such trimers, one could investigate the use 
of this two-level system for quantum computing purposes 
similarly to what was done for quantum qubits using per- 
sistent currents |25j. 
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APPENDIX A: BLOCH SPHERE 

It is convenient to use the Bloch sphere parameteriza- 
tion of the degenerate ground state \g) = (a, /?)* = a\k = 
— 1) + f3\k = 1). Define the unit vector n — (g\<r\g). 
The manifold of degenerate states of the two-level sys- 
tem (undistorted trimer) is represented then by two- 




FIG. 6: In the Bloch sphere parameterization the ground 
states of two-level system correspond to the points of a unit 
sphere i.e., they can be labeled by a unit vector n. For the 
undistorted trimer the xy projection of this unit vector is pro- 
portional to the electric dipole moment while its z component 
is proportional to the magnetic moment of the ground state. 

dimensional sphere swept by the unit vector n. The elec- 
tric dipole moment jl in the corresponding ground state 
is proportional to the projection ri± of ft on the xy plane 
fx = (j, n±. 

The partly-occupied states \k = ±1) of the static 
trimer have an orbital magnetic moment perpendicular 
to the plane of the molecule which couples to a magnetic 
field B. In the degenerate E' subspace, the coupling is 

Hb = -Hoib(TzB z , (Al) 

where the moment /i or b is 2neujA as previously dis- 
cussed. In the Bloch sphere parameterization, the mag- 
netic dipole moment is proportional to the projection of 
n onto the z axis, and is equal to /i rb"z- 

The value of the moment ^ or b for Na3 is about 0.3 Bohr 
magnetons, if we ignore the barrier to free pseudorota- 
tion. The electrical moment ero/2 is 0.98 eA or 4.7D (1 
Debye unit = 0.208 eA.) in our naive model. An unpub- 
lished density functional calculation [26| gives 2. 3D. It is 
not surprising that the actual result is reduced, because 
the true molecular orbitals surely have a significant pa 
component missing from our model. Also, the on-site 
Coulomb repulsion causes a suppression of the dipole as 
is shown in the next appendix. 

APPENDIX B: EFFECT OF A HUBBARD TERM 

The Coulomb interaction between electrons is approx- 
imated by a Hubbard term of strength U, 

H lnt =U ^2 ( c !iT c ™T + c rV«i) ( B1 ) 

n=0,l,2 

The Hilbert space has dimension 20. There are 8 basis 
states with no vacant sites and 12 with one vacant site. 
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with angular momentum k = 1, three with k = — 1, and 
two with k = that have higher energy. Parity symmetry 
is not broken by the Hubbard term. The ground state is a 
k = ±1 doublet, as it was for the non-interacting model. 
Figure [7] shows the dependence of the quantum dipole 
on the dimensionless Hubbard parameter n = U/to. A 
factor of two suppression occurs when U is of the same 
order as to- For k< 1, the doubly degenerate ground 
state energy is E w — 3to + (13/3)[/. The dipole operator 
has an off-diagonal matrix element (— l\fl\ + 1) between 
k = — 1 and k — +1 states equal to — — iy) with 



4 

1 - X K 
9 



for k < 1. 



(B2) 



FIG. 7: The quantum dipole in units ero/2 as a function of 
the dimensionless Hubbard strength re = U/to. The intercept 
coincides with the non-interacting Hiickel model, while for 
K S> 1 it falls as re -3 . 



This agrees with the non- interacting result, Eq. J^J. For 
k ^> 1, the doubly degenerate ground state has energy 
E ps 3£7 — "Hq/U . The magnitude of the off-diagonal 
matrix element falls rapidly as shown in Fig. going 
asymptotically as 



Total spin, spin ^-projection, and orbital (electronic) an- 
gular momentum are conserved. In the eight-dimensional 
5=1/2 and 5 2 = 1/2 subspace, there are three states 



^o, for « > 1- 



(B3) 
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